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We study the relaxation of a spin I that is weakly coupled to a quantum mechanical environment. 
Starting from the microscopic description, we derive a system of coupled relaxation equations within 
the adiabatic approximation. These are valid for arbitrary / and also for a general stationary 
non-equilibrium state of the environment. In the case of equilibrium, the stationary solution of 
the equations becomes the correct Boltzmannian equilibrium distribution for given spin I. The 
relaxation towards the stationary solution is characterized by a set of relaxation times, the longest 
of which can be shorter, by a factor of up to 27, than the relaxation time in the corresponding Bloch 
equations calculated in the standard perturbative way. 



Nuclear magnetic resonance is a well-established method 
for testing electronic properties in solids jlj]. In recent 
years, it became possible to apply this technique not 
only in three dimensions, but also to a two-dimensional 
electron system, the quantum Hall ferromagnet that is 
realized in semiconductor heterostructures in a strong 
magnetic field. The experimental work lead to the unex- 
pected conclusion that a new kind of low-energy states, 
Skyrmions, can be formed and can determine the nuclear 
relaxation processes in these systems when one Landau 
sub-level of one spin direction is filled P~pd| . 

As a theoretical description of spin relaxation, Bloch's 
equations have been successfully used for about fifty years 
now. While these phenomenological equations are appli- 
cable in a wide range of cases, their microscopic deriva- 
tion reveals two main restrictions. First, as was already 
discussed in the original work p^ |, the derivation be- 
comes strictly valid, if either the spin is I = 1/2 or the 
temperature of the bath is large compared to the reso- 
nance frequency. But the spin in the system under study 
can be I = 3/2 (for 69 Ga, 71 Ga, and 75 As; cf. Ref. §), 
or higher in the case of magnetic impurities. Further, 
the progress in the experimental techniques now lets a 
regime of temperatures and magnetic fields come into 
reach, in which the temperature of the bath is of the same 
order as the nuclear resonance energy (nuclear Zeeman 
energy). The second restriction in the derivation of the 
phenomenological equations demands that the environ- 
ment (bath), the quantum mechanical degrees of freedom 
causing the spin relaxation, be in thermodynamic equi- 
librium. But in the case of the quantum Hall ferromag- 
net, the nuclear spins are coupled to a two-dimensional 
electron gas in which the electron-electron interaction 
plays the key role, since all single particle states are de- 
generate into a single Landau level due to the strong 
magnetic field; such a system, dominated by Skyrmion 
states, is not necessarily in equilibrium. Thus, it appears 
worthwhile to reconsider the microscopic derivation of 
phenomenological equations for the spin relaxation in or- 
der to investigate whether there is a significant difference 
between the general case and a case in which the Bloch 
equations are valid. 



It is the purpose of this work to relax the two above re- 
strictions by investigating the general case of an arbitrary 
spin I and also an arbitrary stationary state of the envi- 
ronment responsible for the relaxation. We refrain from 
studying a specific mechanism and consider instead the 
general case of a magnetic moment coupled to a bath of 
other quantum degrees of freedom. This magnetic mo- 
ment can be a nuclear spin, or also a magnetic impurity. 
In the following, we use the terms "nuclear spin" for the 
magnetic moment and "electrons" for the bath - usually 
the latter is called "lattice". Then, the contribution of 
the nuclear spin to the Hamiltonian is 

H = - 7 / ■ (B + B) . (1) 

Here, the magnetic moment, 7/ (where / is the spin), 
couples to an effective magnetic field. It is well known 
that for I > 1/2, there is an additional term caus- 
ing relaxation, the electric quadrupole moment of the 
nucleus coupled to an inhomogeneous external electric 
field. Here, we have omitted this term in H, since the 
model Eq. ([!]) already suffices for the demonstration of 
our method; an inclusion of a quadrupolar coupling is 
straightforward. The effective magnetic field acting on 

the nuclear spin I in Eq. ([l]) contains an operator B 

of the electronic system. One can picture B as being 
proportional to the electrons' spin. Its longitudinal com- 
ponent, B z , modifies the eigenvalues of the nuclear spin 
system, while the transverse components cause transi- 
tions between eigenstates. There is also a fixed part of 
the magnetic field in z-direction, Bq — Boe z , which acts 
as an external field. The coupling beween nuclear spin 
and electrons is supposed to be weak in the sense that 
we can use the adiabatic approximation as discussed be- 
low. We do not need to make any assumptions about the 
electronic subsystem's Hamiltonian or the electronic sub- 
system's state. This Hamiltonian may contain electron- 
electron interactions and the subsystem may be in an ar- 
bitrary stationary state, equilibrium or non-equilibrium. 
We want to derive kinetic equations for the expectation 

value of the spin vector I (I = 1(1 + 1)). To this end, we 
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need to take all components of the density matrix into 
account, not just only the spin vector. Then, as will be 
seen below, it is advantageous to use the spherical tensor 
operators Tlm as a complete basis in the space of oper- 
ators acting on the state of a spin /. The Tlm are the 
irreducible tensor operators in the spherical coordinate 
representation |L3]]. For actual calculations, the follow- 
ing definition using the spherical harmonics Ylm proves 
very helpful: 



f LM := Mil (I-V) l r L Y LM (f) 



(2) 



Tlm is a polynomial in the components of the spin oper- 
ator. It is independent of the auxiliary variable r, since 
on the r.h.s., there are L derivatives acting on a poly- 
nomial of order L (f = f/r denotes the unit vector; we 



use the conventions of Rcf. 13 for the spherical harmon- 
ics Y L m)- For a spin /, the (21 + l) 2 operators T L m 
with L — ■ ■ • 21 and M = —L ■ ■ ■ L form a complete 
system of operators acting in the spinspace. The nor- 
malization, AfiL = 2 L v/47r(27 - L)\/(2I + L + Tfi/Ll, is 

choosen such that Sp j T| M f L > a/'| = 5 L l> 5 m m' , 

and we have f\ M = {-l) M T L ,- M - Specific expressions 
for the operators Tl,m (f° r L = • • • 4) can be found in 
Table IV of Ref. Q'. 

After having established the basic notation, we proceed 
now to describe the derivation of the kinetic equation for 
the average Tlm{^) = (Tlm), to find the stationary solu- 
tion and to study finally the relaxation towards the sta- 
tionary solution. Here and below, the brackets (• • •) stand 
for the state of the combined system of spin and electrons. 
We use the framework of the Keldysh method |15| in or- 
der to derive the kinetic equation. As in our earlier work 
on the electron spin relaxation Jl6| , we shall employ the 
adiabatic approximation; that means that in the equa- 
tion of motion, the effect of the coupling between spin 
and electrons is neglected beyond the first order in the 
spin's eigenenergies and also neglected beyond the second 
order in the relaxation times, cf. Ref. |17|] . The unper- 
turbed motion of the spin is a precession with the fre- 
quency u>o = 7 Bq. Then, we have up to the second order 

of the perturbation theory in the coupling V — —7/ • B: 
(id t - M lu ) T LM (t) = 



T LM (t),V(t) 



f LM (t),V(t) ,V(h 



0{V 3 



(3) 



In the spirit of the adiabatic approximation, we now 
consider the terms of the perturbation series in higher 
than second order as giving rise to an additional (weak) 
time dependence of the spin operators due to the relax- 
ation process, and we decouple the expectation values 
of spin and electron operators. The first order term, 



{[T LM ,I](i)) ■ (B(t)), describes the Knight shift, the shift 
of the nuclear resonance frequency due to the coupling to 
the electrons. Since we want to focus on the relaxation, 
we disregard this corrections of the spin's eigenfrequency 
in the following. Then, we get 



(t d t - Mlj ) T LM {t) 



£ (« 



rn-.tn' — 0.±1 



1 *n 
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Tlm, In 



'}(<)> 



R—m\ — n 



(4) 



In deriving Eq. (||) , we took the unperturbed time depen- 
dence of the operator I(t\) into account, but neglected 
the difference t — t\ in its time dependence due to the 
relaxation. In Eq. (Q), equal-time commutators ([,]) and 
anticommutators ({, }) can be evaluated, and the expec- 
tation values can again be expressed by Tlm^) as we 
shall see below. Thus, it is the use of the spherical tensor 
operators that makes it possible to derive closed coupled 
equations for the relaxation of the spin. Most crucial 
are the averages of the electronic subsystem that enter 
Eq. (|J). Since we do not assume thermodynamic equi- 
librium for the electronic subsystem, we get both, a cor- 
relation function and also a response function which are 
independent and which we denote by C m - m i and R m ;m' , 
respectively. They are given by 



dt\e 



'^''-^({BmW-^ttl)}) (5) 



-im' ujq (t — ti ) 



B m (t),B TO /(ti) ) (6) 



Our convention for the vector components of / and B 
is: J±i = 4 ± il y , Iq = I z and B ±1 = (B x ± iB y )/2, 
Bo = B z . 

Commutators and anticommutators in Eq. (Q) can be cal- 
culated either directly from the definition Eq. (|^) , or with 
the aid of the theory of spherical tensors (cf. Ref. |^8|). 
The commutators express the behavior of the Tlm under 
rotations: 



TlmJz — —M Tj 



LM 



(7a) 



Tlm, I± 



= -y/L(L + l)-M(M±l) T l ,m±i . (7b) 



The anticommutators can also be expanded in Tlm] spe- 
cializing the general result containing Racah coefficients 
(see Eq. (17) of ref. 0) to our case, we get for m = 0, ±1: 



^T L M,im^ — 0>L+l-M TL+l,M+m + b^lj f L -l,M+r 



(8) 
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with the coefficients 

«L;M = &i;M = W-M* C L , (9) 

«S = &EL = T V(L + l±M)(L±M) c L , (10) 
where 

/ (21 +1)2- £2 

Ci = V (2L) 2 -1 • (H ) 

Inserting relations @) and (||) into Eq. (Q) solves our 
task of deriving a closed set of relaxation equations for 
arbitrary nuclear spin I. The equations are linear in the 
expectation values of the spherical tensor operators Tlm- 
The properties of the electronic system enter the equa- 
tions parametrically in the form of correlation functions 
and response functions. Since the spin vector is given by 

the Tim, a relaxation equation for (I) can be extracted 
from the system of equations (Q). In the case of I = 1/2, 
this is particularly simple, since Tlm = for L > 1 (and 
Too = 1/V2 is constant) and we recover the Bloch equa- 
tions, see below. 

In this general form, the equations (^) are still not very 
transparent. Therefore, we now make additional assump- 
tions regarding the correlation and response functions en- 
tering Eq. (^) . The term I z B z in the perturbation V 
(corresponding to m, ml = in Eq. (^)) just changes the 
spin's resonance frequency luq. Since we already omitted 
the Knight shift, the first order correction term to this 
frequency, we now also omit consistently the second order 
contributions resulting from m = or m! = in Eq. (|J). 
Next, the terms with m — m' in Eq. (Q) are neglected 
too. In the electron system, they would correspond to a 
total change of the z-component of the spin by 2; if the 
electrons' state is a strict eigenstate to z-component of 
the spin, such expectation values vanish. Under these 
assumptions, both correlation and response functions, 
Cm\m' and Rm-m' , now become "diagonal" in the index 

^5 Cm;m' ^m; — m'Cm and Rm,\m' — ^m;— m'-^m* Now, 

from the general structure of Eq. (|J), it is obvious that 
the equations do not couple Tlm(^) for different M. So 
finally, we arrive at our general result, valid for arbitrary 
spin I and for a non-equilibrium state of the bath 

(d t + iMuuo) T LM {t) = 

[i(T+l)-M 2 ] T LM (t) 

{l^(L + 1) 2 - M 2 c L+ i T L+hM {t) 

- (L + 1)^L 2 - M 2 c L T L ^ liM (t)} . (12) 

We have defined two different times by 1/t + := 2(C\ + 
C-i) and 1/t_ := 2(i?i— -R-i). These times are indepen- 
dent, as long as we assume a non-equilibrium state of the 



electronic system. We have r+ > and the ratio t + /t_ 
can be expressed as t + /t_ = (1 — q)/(l + q) with positive 
q. Thus, |t+/t_| < 1. If the electrons are in equilibrium 
at temperature T, the fluctuation-dissipation theorem 
results in q = exp(— uq/T). 

We want to stress that these relaxation equations are 
valid (under the assumptions stated above), no matter 
which specific relaxation mechanism one wants to con- 
sider. The special mechanism enters the equations in the 
form of two time scales, t + and t_. In the case of ther- 
modynamic equilibrium in the bath, their ratio, t + /t_, 
is fixed by the temperature, and only the time r + - which 
also depends on temperature - is specific for the relax- 
ation mechanism. 

Relaxation equations generally serve two purposes. First, 
they determine a stationary state. Second, they describe 
the relaxation towards this state. We now want to discuss 
both of these points. 
Stationary solution 

The study of the stationary solution of Eq. ([[2]) serves 
as an important test of our procedure, since, in the case 
of equilibrium in the bath, the result is obvious. With 
the ansatz Tlm(^) = ^m-,o T^, one derives a recursion 
relation 

CL+i T L+1 = -— T L + c L Tz_i , L = l---27 (13) 

T+ 

where T = (21 + 1)" 1/2 , T 2J+1 = 0. The solution is 

h = \l{I+l){l-q)Mi X /Mix (14) 

where Mij are elements of the matrix Ai = M\ ■ 
M 2 ■ ■ ■ M11 with 

m, = ( 1+ ; v-f^). da. 

Evaluation of this matrix product yields - for a general 
non-equilibrium state of the bath - the following station- 
ary distribution for the z-component of the spin vector 

I z = ]T m q- m I ]T q- m . (16) 
m——I m——I 

If we assume now the electronic system is in equilibrium 
at temperature T, we find that I z is given by the Brillouin 
function, the well-known correct equilibrium distribution 
for a spin I. 
Relaxation 

The general solution of the coupled relaxation equations 
Eq. ( |l2] ) is a superposition of exponentially decaying 
terms. We are mostly interested in the relaxation of the 
spin components I z oc T\q and (I x + il y ) oc Tu, which 
are described by the part M = and M = 1 of the full 
system of equations. For M = and M = 1, there are 
21 different relaxation times, for M = 2 there are 27—1 
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times; in total we have 1(21 + 3) times. Here, each case 
of I needs to be discussed separately. 
We start with spin 1/2. In this case, the En. ([lj) gives 
immediately the standard Bloch equations M : 



(a f + itjo) /+(t) 



(17) 




FIG. 1. Spin I = 3/2: Solution of the relaxation equations 
for M = with t+/t_ = 0.8; 7 Z oc Xio (arbitrary units used). 

More interesting, because of its experimental relevance 
||, is the case of spin 3/2. The relaxation equations for 
M = couple the expectation values of three operators: 
Tio, T20, and T^o- These equations are easily solved. In 
Fig. [IJ we show the flow of the expectation values to- 
wards the stationary state at the origin in the case of 
t + /t_ = 0.8. The inverse time scales, A, which deter- 
mine the relaxation are shown for both the longitudinal 
(Aj,, M — 0) and the transverse equations (At, M = 1) in 
units of 1 /t + as functions of the temperature in the case 
of equilibrium in the bath in Fig. |[ Between the high- 
temperature limit T » luq and the low-temperature 
limit T << luq, both the largest longitudinal and the 
largest transverse relaxation time decrease by a factor of 
three as compared to the relaxation time in Bloch ap- 
proximation calculated in the usual perturbative way. 
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FIG. 2. Spin I = 3/2: Inverse relaxation times in the lon- 
gitudinal (solid line) and transverse (dashed line) relaxation 
equations in units of l/r+. 

For general /, we determine both the longitudinal and 
the transverse relaxation times in the following limit- 
ing cases. For T >> u (t+/t_ ~ 0), the longitudi- 
nal relaxation time is t + and the transverse 2t + as in 
the Bloch equations. In the opposite case, T « uj 



(t+/t_ — > 1), all the operators Tlm become equally im- 
portant in Eq. (|l2|). Then, the 21 relaxation times are 
as follows (n = 0, 1, • • • 21 — 1): In the longitudinal case, 
M = 0, we get r + /[2(n + 1)1 — n(n + 1)] and the longest 
relaxation time (n = 0) is always twofold degenerate; in 
the transverse case, M = 1, we get r+/[(2n + 1)1 — n 2 ] 
and here, the longest relaxation time (n — 0) is always 
non-degenerate. These explicite expressions are a con- 
jecture based on an explicite evaluation of the relaxation 
equations for spin I up to 7/2. 

Starting from the microscopic description, we have de- 
rived relaxation equations valid for an arbitrary spin / 
coupled to a quantum mechanical environment in a non- 
equilibrium state. The solution of these coupled equa- 
tions shows that, compared to the Bloch equations, there 
is an additional temperature dependence in the relax- 
ation time which can decrease the relaxation time by a 
factor of up to 21. 
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